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The dispersion of index-guiding microstructured polymer optical fibers is calculated for second-harmonic 
generation. The quadratic nonlinearity is assumed to come from poling of the polymer, which in this study 
is chosen to be the cyclic olefin copolymer Topas. We found a very large phase mismatch between the pump 
and the second-harmonic waves. Therefore the potential for cascaded quadratic second-harmonic generation is 
investigated in particular for soliton compression of fs pulses. We found that excitation of temporal solitons 
from cascaded quadratic nonlinearities requires an effective quadratic nonlinearity of 5 pm/V or more. This 
might be reduced if a polymer with a low Kerr nonlinear refractive index is used. We also found that the 
group-velocity mismatch could be minimized if the design parameters of the microstructured fiber are chosen 
so the relative hole size is large and the hole pitch is on the order of the pump wavelength. Almost all 
design-parameter combinations resulted in cascaded effects in the stationary regime, where efficient and clean 
soliton compression can be found. We therefore did not sec any benefit from choosing a fiber design where 
the group-velocity mismatch was minimized. Instead numerical simulations showed excellent compression of 
A = 800 nm 120 fs pulses with nj pulse energy to few-cycle duration using a standard endlessly single-mode 
design with a relative hole size of 0.4. © 2008 Optical Society of America 
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1. Introduction 

Microstructured optical fibers allows for an elaborate 
dispersion control. Usually they are exploited to gain 
control over the group- velocity dispersion (GVD) [1] (see 
[2] for a recent literature overview over the dispersion 
properties of microstructured optical fibers). Also the 
group-velocity mismatch (GVM) for three-wave mixing 
processes can be completely removed [3], which is im- 
portant in order to realize second-harmonic generation 
(SHG) of ultra-short fs pulses. However, the prize for 
having a small GVM is that the phase mismatch be- 
comes very large [3], and efficient SHG would therefore 
rely on quasi-phase matching (QPM) techniques. 

The large phase-mismatch A/3 can instead be ex- 



ploited for cascaded x 



(2) 



x(2) 



SHG processes [4] 



The second harmonic (SH) is within a coherence length 
27r/|A/3| generated and then back-converted to the fun- 
damental wave (FW). This cascaded process generates 
a nonlinear phase shift <I>nl on the FW [5], which can 
become several units of tt by propagating in just a few 
em's of nonlinear material, and is conceptually equiva- 
lent to the nonlinear phase shift observed through self- 
phase modulation (SPM) with cubic nonlinearities [6]. 

The large nonlinear phase shift can be exploited to 
compress ultra-short fs pulses [7]. This is particularly 
advantageous when the compression occurs inside the 
nonlinear material due to the soliton effect [8,9]; the 
cascaded SHG process can generate a negative $nl re- 
sulting in a negatively chirped FW, and if the mate- 
rial has normal GVD a temporal soliton is generated. 
By launching a higher-order soliton, the input pulse can 
be compressed by exploiting the initial pulse narrowing. 



This cascaded quadratic soliton compressor (CQSC) can 
in principle compress the FW pulse to single-cycle dura- 
tion, and is ultimately limited by higher-order dispersion 
and competing Kerr nonlinear effects [10]. 

One requirement for efficient compression is that GVM 
effects are not too strong: they tend to distort the com- 
pressed pulses through a Raman-like effect [9-12]. In 
fact, when GVM dominates over the cascaded effects 
from the phase mismatch the compression is nonstation- 
ary, resulting in inefRcient compression and distorted 
pulses. Therefore the possibility offered by microstruc- 
tured optical fibers to control GVM is very intriguing 
for the CQSC, because the compression can become sta- 
tionary, which implies efficient compression and clean 
pulses. The fiber geometry can also help overcoming the 
problem of inhomogeneous compression in the transverse 
direction of the beam found in a bulk geometry [13] 

We have done a preliminary investigation of the po- 
tential in using silica microstructured optical fibers for 
CQSC [2,14], where the quadratic nonlinearity was as- 
sumed to come from thermal poling of the silica fiber 
[15]. We surprisingly found that zero GVM was not as 
such an advantage, since the compressed pulses became 
very distorted. Another main result was that a very large 
quadratic nonlinearity deff '--^ 3 — 5 pm/V was needed in 
order to generate suitably large $nl [2] . This is because 
the chosen fiber designs had very a large phase mismatch, 
and <I>NL oc Realistically, thermal poling of sil- 

ica fibers can generate dcff ~ 0.5 pm/V [16]. Therefore 



^No comparison should otherwise be made between the bulk 
and the fiber CQSC; the bulk version works with high-energy fs 
pulses with /.ij-mj energies, while the fiber version works with low- 
energy fs pulses with pj-nj energies. 
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we suggested to lower the phase mismatch using QPM 
techniques, and in this way a much lower deft ~ 1 pm/V 
was sufficient [14]. However, periodic thermal poling of 
silica fibers have until now proved inefficient, and has 
resulted in a much lower doff ^ 0.01 pm/V than ex- 
pected [17]. We must therefore conclude that periodi- 
cally poled silica (microstructured) fibers cannot gener- 
ate large enough <i>NL to be interesting for CQSC. 

A large doff would compensate for a large A/3. In this 
paper, we therefore turn our attention to microstruc- 
tured polymer optical fibers (mPOFs). Poling of poly- 
mers can generate extremely large quadratic nonlineari- 
ties (ranging from 1 pm/V to lOO's of pm/V) [18], while 
poling still has to be shown in a fiber context. Alter- 
natively, the low mPOF drawing temperature (few hun- 
dred degrees C) implies that nanomaterials with large 
quadratic nonlinear responses (such as nanotubes [19]) 
can be drawn into the fiber core. Also this solution 
promises a strong doff. The potential for using mPOFs 
for quadratic nonlinear optics is therefore large. 

Here we calculate the dispersion properties of an 
index-guiding mPOF with three rings of air-holes in the 
cladding. The polymer material is chosen to be the cyclic 
olefin copolymer Topas due to its broad transparency 
window [20]. We point out that a problem with using 
polymer as fiber material is that unlike silica the ma- 
terial Kerr nonlincarity is very large, and it generates 
an SPM-induced positive nonlinear phase shift which 
has to be overcome by the negative cascaded nonlin- 
ear phase shift. We will show that for the considered 
polymer Topas doff > 5 pm/V is needed to do so (com- 
pared to doff ~ 3 — 5 pm/V in silica [2, 14]). Such a value 
could be achieved with polymers as fiber material. We 
also show that the fiber designs with a dramatically re- 
duced GVM are multi-moded in the SH, and thus one 
risks to have cascaded nonlinear conversion to several SH 
modes. Quite surprisingly we find that the compression 
is always in the stationary regime, irrespective of the 
choice of fiber-design parameters. This is very positive 
since efficient and clean compression can be obtained. 
On the other hand, there is no longer any motivation for 
reducing GVM in order to enter the stationary regime. 
Thus, we conclude that there are very few benefits of 
reducing GVM through fiber design. This statement is 
underlined by performing numerical simulations of the 
propagation equations for an mPOF having doff = 10 
pm/V: an endlessly single-moded mPOF design dom- 
inated mainly by material dispersion turns out to give 
clean and efficient compression of fs pulses, and obtaining 
compressed pulses with durations close to single-cycle 
duration is possible. So if the potentially large doff of 
polymer material can be exploited, then excellent pulse 
compression can be obtained in mPOFs. 



2. Nonlocal theory 

The scope of this theoretical part is to point out the im- 
portant parameters when choosing the proper fiber de- 



sign. The main hypothesis is that an mPOF can be used 
as a CQSC. The challenges with creating a quadratic 
nonlincarity in the fiber aside, there are other obsta- 
cles that must be faced. Understanding these issues can 
be greatly enhanced by realizing that in the cascading 
regime, the coupled FW and SH propagation equations 
reduce to a nonlinear Schrodinger equation (NLSE) for 
the FW [6]. The action of the cascaded SHG can be 
modeled as a Kerr-likc nonlincarity with an equivalent 
nonlinear refractive index tl^^q cx d^g/A/3. Importantly, 
when A/3 = /32 — 2/3i > (which is usually the case) n^i^Q 
is self-defocusing of nature: it generates a negative non- 
linear phase shift $nl < 0. The material Kerr nonlinear 
refractive indejH ?^Korr instead usually self-focusing, 
and therefore counteracts the effects of the cascaded non- 
linearities. Achieving n^^Q > n^g^.^. is crucial to obtain a 
large negative nonlinear phase shift. An important point 
here is that because the cascaded SHG induces a self- 
defocusing Kcrr-like nonlincarity, temporal solitons ex- 
ists in presence of normal GVD. In contrast, the usual 
cubic temporal solitons observed for instance in telecom 
fibers require anomalous GVD due to the self-focusing 
Kerr nonlincarity of silica [21]. 

GVM is also playing a decisive role. As we recently 
showed, the pulse compression is clean and efficient in 
the stationary regime, where the phase mismatch dom- 
inates over GVM effects [10, 12]. These results were ob- 
tained by deriving a more general NLSE, where disper- 
sion including GVM imposes a Raman-like nonlocal (de- 
layed) temporal cubic nonlincarity. Contrary to earlier 
reports [7, 11] we showed that the regime where effi- 
cient compression takes place is independent on the input 
pulse duration. 

In the remainder of this section we briefly show how 
the nonlocal NLSE for the FW is derived, and discuss the 
consequences for soliton pulse compression. The equa- 
tion is derived in dimcnsionlcss form based on the SHG 
propagation equations derived in App. [X] 

The nonlocal theory was derived in Refs. [10, 12] (for a 
general review on nonlocal effects, see [22]). Essentially 
in the cascading limit (large phase mismatch), the SH 
becomes slaved to the FW. The normalized FW Ui can 
then on dimensionless form be modeled by the following 
nonlocal nonlinear Schrodinger equation (NLSE) [12] 



Ui+N^,,,Ui\U,\^ (1) 



2 9r2 

pOO 



where Kerr cross-phase modulation (XPM) effects have 
been neglected, and for simplicity self-steepening and 
higher-order dispersion are not considered. On the left- 
hand side an ordinary NLSE appears with an SPM-term 
from self- focusing material Kerr nonlinearities. On the 
right-hand side the effects of the cascaded SHG appear: 



^The Kerr nonlinear refractive index is usually denoted n2, but 
we reserve this subscript to the SH. 
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it is also SPM-like in nature, but it turns out to be con- 
trolled by a temporal nonlocal response. The dimension- 
less temporal nonlocal response function R± (r) appears 
in two distinct ways: in the stationary regime it is given 



by and \R+\ oc e 



is localized having a char- 



acteristic width Tf, of a few fs. This happens when the 
phase mismatch is dominating over GVM effects, or more 
precisely when A/3 > A/3sr, with 



and I?2 = ^Stt/Jj^-'/Aj is the fiber GVD parameter. In 
the stationary regime clean and efficient compression can 
be obtained [9,10,12,23]. Instead in the nonstationary 
regime (A/3 < Ap^r) GVM effects dominate, and the 
compression becomes distorted and inefficient. In this 
case the nonlocal response function is given by i?_ , and 
cx sin(|r|/rb) is oscillatory and never decays. 
The fiber designs presented here all turn out to 
be in the stationary regime. In this case, and when 
the nonlocal response function can be assumed quasi- 
instantaneous, Eq. (HI) can be written as [9,11,12] 



(2) 
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The soliton order for the Kerr fiber nonlinearity is [24] 



Kerr 



^Kcr 



(4) 



where Ld,i = Ti^i/I/^P^I is the characteristic GVD length 
of the FW, and Lkbti = (7Kerr^in)~"^ is the characteristic 
Kerr nonHnear length. Pin = eoneS,icai£^^/2 is the input 
peak power, related to the peak input electric field 
The Kerr nonlinear coefficient is 



7Ke 



UJl 



'Kerr 



cA. 



eff,l 



(5) 



where ^eff,i is the effective Kerr FW mode overlap area 



_ (/rfx|fi(x)|^)^ 



(6) 



Note that Acff.i = /j^^ from Eq. and that as ex- 

plained in App. |X]the subscript 'P' indicates that inte- 
gration must be done only over the part of the fiber 
mode sitting in the polymer part of the mPOF. We 
also encounter the well-known Kerr nonlinear refrac- 
tive index f^Kcrr• ^^^^ refractive index change due to 
the Kerr self-focusing effect is defined as An = n^.^^^!, 
where / is the intensity of the beam. Although it is not 
known for the Topas polymer, we estimate (based on 
similar polymer types) that it is in the range n^.^^.^ ~ 



10 — 15 • 10^^" m^/W, which is an order of magnitude 
larger than for fused silica. 

In a similar way we have introduced a soliton order 
from the cascaded SHG process 
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where Xshg = (TSHG^ln) ^- Here the nonlinear coeffi- 
cient is defined as for the Kerr case 
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only it contains the SHG overlap area Ashg 
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/pdx[F;(x)]2F2(x)|^ 
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where Oo are the fiber mode areas (IT9l). The "effective" 



nonlinear refractive index from the cascaded process [5] 
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Note the negative sign in front of ngjjQ; since the cases 
we will discuss always have A/3 > the cascaded non- 
linearity is therefore self-defocusing. Thus, in order to 
generate a negative nonlinear phase shift we must have 
7SHG > 7Kcir- We can express this through an effective 
soliton order 



NIhG - ^Kcrr = ^in^D,! (7SHG " 7Kcrr) (H) 



Using this in Eq. ^ we then get an NLSE with a self- 
defocusing SPM term: 
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Ui-N^sUi\Ui 



= -lNl^aSaTR,SHG\Ul\^-^ (12) 

Self-defocusing solitons with strength N^ff can then be 

(2) 

excited if the FW GVD is normal, i.e. if /3} ' > 0, which 
has been taken into account in Eq. (H^). 

So the LHS of Eq. (HH) is now a self-defocusing NLSE 
supporting solitons if TVcfr > 1. The RHS contains the 
Raman-like perturbation due to GVM. It becomes im- 
portant for large SHG soliton orders iVsHG [10] and 
when the compressed pulse duration is on the order of 
7fl,SHG = T'iJ.SHGTin, where the characteristic dimcn- 
sionless time scale of the Raman-like perturbation is 



2\d, 



TR,SliG 



Am. 



(13) 



Typically Tr^shq is on the order of 1-5 fs, but can be- 
come very large in the nonstationary regime [10]. Finally, 
(2) 

Sa = sgn((ii2/32 ): the Raman-like effect will depending 
on this sign give either a red-shifting {sa = +1) or blue- 
shifting {sa = —1) of the FW spectrum. In the cases we 
present here Sa = —1. 
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Fig. 1. The microstructured fiber design considered here. 
3. Numerical results 

In this section the transverse fiber modes and their dis- 
persion properties are calculated for various mPOF de- 
signs (varying the air-hole diameter d and the hole pitch 
A). Then we show compression examples for selected 
fiber designs. These numerical simulations were done us- 
ing the full coupled propagation equations (|22|) , which in- 
clude Kerr XPM effects, steepening terms, and are valid 
down to single-cycle resolution using the slowly-evolving 
wave approximation [23,25,26]. 

A. The microstructured fiber 

We consider an index-guiding mPOF with a triangular 
air-hole pattern in the cladding, see Fig. [TJ The mPOF 
design parameters are the pitch A between the air holes 
and the relative air hole size D = d/A, where d is the 
physical air-hole diameter. We assume that the core has 
a quadratic nonlinearity from thermal poling of the poly- 
mer after the fiber drawing, or from including nonlinear 
nanomaterial in the fiber preform. 

We consider an mPOF made from the cyclic olefin 
copolymer Topas [20]. The advantage of using Topas 
compared to, e.g., poly-(methyl methacrylate) (PMMA), 
is that Topas has a larger transparency window: it is 
transparent from 300-1700 nm, interrupted by two ab- 
sorption peaks around 1200 nm and 1450 nm [20], see 
App. [Clfor more details. 

B. Fiber design, dispersion and compression 

We will here search for fiber design parameters, where 
the requirement is an optimal dispersion profile for cas- 
caded quadratic soliton compression. Details about the 
fiber mode calculations and the dispersion calculations 
arc found in App. |B1 

Two types of designs are discussed: one is a realistic 
mPOF with a large air-hole pitch of around A = 7 /xm 
that can easily be drawn, while another takes A small 
(comparable to the FW wavelength) as to significantly 
alter the dispersion parameters. Such an mPOF would 
probably need to be done by tapering an mPOF with a 
larger pitch. In the calculations the Sellmeier equation 
from App. [O with T = 25°C was used. 

It is fairly standard to draw mPOFs with a air-hole 



Fig. 2. The GVM parameter di2 and phase mismatch 
A(3 for an mPOF with A = 7.0 /.im. In (a) the fines for 
single-mode operation are indicated for the FW and the 
SH modes. The FW zero dispersion point is also indi- 
cated; below this line the FW GVD is normal. Finally 
the 'X' indicates the chosen fiber design for single-mode 
operation at Ai = 800 nm. 



pith of A > 5 /im. With this as a starting point we show 
in Fig. [2] the GVM and phase mismatch for a Topas 
mPOF with A = 7.0 iim. We notice that the choice of 
D does not affect the dispersion parameters much. This 
is because the size of the waveguide [the core diame- 
ter is dcore = A(2 — D)] is significantly larger than the 
wavelength of the guided light, and so material disper- 
sion is dominating (see also Fig.[3|). The microstructured 
cladding therefore alters the dispersion only very little, 
and we cannot get tailored GVM or GVD properties. 

Since Topas has a primary transparency window be- 
tween 290-1210 nm, we focus now on SHG with Ai = 800 
nm. We would like the fiber to be single mode both 
at the FW and at the SH wavelength, because then 
phase matching can only occur between these modes 
and we do not have problems with the typical higher- 
order mode interaction usually observed in wave-guided 
SHG. The criterion for single- mode operation in the 
particular microstructured fiber, we investigate, is [27] 
A/A = 2.8(1? - 0.406)° *^^. This criterion gives the lines 
shown in Fig. ^a.) for the FW (where A = Ai) and the 
SH (where A = A2 = Ai/2). We therefore choose a design 
with D = 0.4, which is indicated with an 'X'; this design 
is actually endlessly single- moded [28]. 

For the chosen design, we show in Fig. [3] the disper- 
sion parameters as function of wavelength. In (a) the 
phase mismatch A/? is shown. It is always above the 
boundary to the stationary regime as given by Eq. (U) 
A/Jsr = di2/2(32^\ Thus, the chosen fiber design for all 
the shown wavelengths is in the stationary regime for 
soliton compression, which is very good news for the 
possibility of clean and efficient compression. This might 
seem surprising given the large GVM seem in Fig. [3Kb), 
but it is a consequence of having a very large phase mis- 
match and the very large GVD at the SH wavelength, 
see Fig. [JJc). However, the large phase mismatch unfor- 
tunately also implies a small cascaded nonlinear param- 
eter 7SHG- In order to achieve 7shg > 7Korr, as required 
for solitons to exist, the effective quadratic nonlinearity 
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Fig. 4. Numerical simulation of soliton compression in an 
mPOF with the design parameters of Fig.[3]and doff = 10 
pm/V. Top left plot shows the FW compressing to 15 
fs FWHM after 30 mm of propagation (at the dashed 
line) starting from 120 fs. Top right shows the spectral 
broadening and development of SPM-induccd sidebands 
of the FW. The pulse parameters were Ncs ~ 3, Pin = 
4.6 kW, and the pulse energy was 0.6 nJ. The middle row 
shows a simulation for a higher soliton order (N^g = 9, 
Pin = 41.8 kW, and a pulse energy of 5.7 nJ), which 
compresses to 4.4 fs after 6 mm of propagation. A cut at 
2 = 6 mm is shown in the bottom plot. Up to 10th order 
dispersion was included {irid = 10). 2^^ temporal points 
and > 15 z-steps per coherence length were used. 



Fig. 3. The dispersion for A = 7.0 /im and D = d/ A ^ 
0.4. (a) the phase mismatch Af3 = P2 — 2/3i and (b) 
GVM di2 = /3f ^ - vs. Ai. (c) the GVD parameter 
V = — 27rc/3'^^'/A^ vs. A. The material dispersion is also 
shown, (d) shows the effective nonlinearity required to 
achieve 7shg > 7Korr for n^^^.^. = 15 • 10"^° m^/W. 



doff must be quite high, around 6 pm/V for Ai = 800 nm 
as shown in Fig. ^d). 

In Fig. |4] we show numerical simulations of soliton 
compression using the chosen design. For the Kerr non- 
linear refractive index of the polymer Topas we use 
'^Kcrr = 15 ' 10~^° m^/W, which is a large but realis- 
tic value. In order to outbalance this Kerr self-focusing 
nonlinearity the quadratic nonlinearity must be doff > 6 
pm/V, see Fig. [3l[d), and doff = 10.0 pm/V was chosen. 
The compression is shown for both a low and a high soli- 



ton order, and both are very clean for several reasons: 
firstly, the compression occurs in the stationary regime, 
and secondly the Raman-like perturbation is very small, 
7j?,,SHG = 1-3 fs. Lastly, looking at the FW frequency 
spectra there are no dispersive waves. These tend to in- 
duce trailing oscillations on the compressed pulse [10], 
but we checked that for this fiber design they are not 
phase matched in the transparent region of Topas. We 
should also mention that wc have neglected any cubic 
Raman effects of the material due to the lack of material 
knowledge in this respect, and if significant Raman ef- 
fects are present as for silica fibers then this would have 
an impact on the results. 

So a standard mPOF, where the fiber simply serves 
as a waveguide but otherwise has little influence on the 
dispersion, can give decent compression (provided doff is 
high enough). What benefits can we get with a reduced 
GVM? To get a detailed understanding of this, let us in- 
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Fig. 5. As Fig. m but with A = 0.65 yum. A zero GVM 
curve is indicated with black as well as the FW zero dis- 
persion point (ZDP). The two 'X'es indicates the chosen 
fiber designs at Ai = 1040 nm. 



vestigate a possible fiber design having very small GVM. 

In order to significantly alter the dispersion, the 
mPOF must have a lower hole pitch. As an extreme ex- 
ample Fig. [S] shows the calculated fiber properties for 
A = 0.65 /xm the GVM can be very small or even 
zero along the black curve. Unfortunately for Ai = 800 
nm this low-GVM area is in an anomalous dispersion 
region, where cascaded solitons do not exist, so let us 
instead work at Ai = 1040 nm, which would be relevant 
for compressing the output of Yb:doped fiber lasers. 

The dispersion for Ai = 1040 nm of this mPOF 
with A = 0.65 /.tm is shown in Fig. [6l^a) as function 
of D. Clearly the dispersion changes dramatically as 
D is changed. The GVM goes from a normal value 
{di2 ~ —200 fs/mm) for low D to zero at high D. 
Nonetheless, we remain in the stationary regime even for 
low D-values since A/3 > A/3sr always. This is because 
both A/? and the SH GVD 2?2 are large. Notice also that 
the FW third order dispersion f3\ changes dramatically: 
it is very large for high D and close to zero for low D. 
Finally, Fig. [6l^b) shows the critical dcs required in or- 
der to overcome the material Kerr self-focusing effects. 
It is around 2-3 times higher than in Fig. [3l which is a 
consequence of a large A/3. 

Zero GVM can therefore be achieved choosing D — 
0.95 [at the cross marked (1) in Fig. [5]. This choice leaves 
the SH multi-moded, and is therefore not suitable for 
compression: there would be different modes with differ- 
ent degrees of phase mismatch interacting with the FW, 
and gaining control over the compression would be dif- 
ficult. If we want the SH to be single moded, we should 
choose a design below D = 0.6 (see Fig. [5]), and the de- 
sign marked with 'X' (2) in Fig. [5] has the advantage of 
being endlessly single moded. In this regime the GVM is 
not too different from material dispersion. 

So the penalty of choosing a small-pitched fiber in 
order to reduce the GVM seems to be too high. The 
main motivation behind controlling GVM through the 
fiber microstructure was to be in the stationary regime. 




^This value is unrealistic but serves the purpose of discussing 
the conditions under which zero or low GVM can be achieved. 



0.6 0.1 

D=d/A 



Fig. 6. (a) The dispersion as function of the relative hole 
size D as calculated for Ai = 1040 nm and A = 0.65 /J,m. 
(b) The the effective nonlinearity required to achieve 
7SHG > 7Korr assumiug = 15 • 10"^" m^/W. 



However, almost all mPOF designs we have investi- 
gated have been in the stationary regime thanks to the 
large phase mismatch. Another motivation was to re- 
duce the characteristic time of the Raman-like response 
T'fl.SHG = M12I/2A/3 as to get cleaner compressed pulses, 
but T/f^sHG is low even when GVM is not reduced again 
thanks to the large phase mismatch. What is more, 
choosing a zero GVM design surprisingly has been shown 
not to benefit compression: In our previous study [14] 
we performed simulations with just the lowest SH trans- 
verse mode (neglecting the higher-order modes), and the 
fiber designs with zero GVM were giving quite distorted 
compressed pulses. This can be partially explained by 
the fact that higher-order dispersion is significant when 
the relative hole size is large; observe for instance how 
large the third-order dispersion is in Fig. [6] for D close to 
unity: this value is several order of magnitude stronger 
than the material third-order dispersion. 
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Fig. 7. The dispersion calculated for D = 0.4 nm and 
A = 1.0 fiia. The lower plot shows the results of a nu- 
merical simulation, where a Ai = 800 nm Tin = 120 fs 
FWHM input pulse is compressed to 5.5 fs. The pulse 
cut is shown after z ~ 7.0 mm of propagation. The fiber 
had dcff = 10 pm/V as in Fig. 31 The input pulse had 
iVoff = 9, Pin = 3.1 kW, and pulse energy 0.42 nJ. 

An advantage of small-pitched fibers is instead that 
nonlinear effects occur with a much smaller power. As 
an example, choosing A = 1.0 /im and D = 0.4 gives 
an endlessly single-mode mPOF, where GVM and GVD 
are similar to the material dispersion (see the calculated 
dispersion parameters in Fig.[7|). This situation is there- 
fore similar to the large-pitch case discussed previously 
(where A = 7.0 However, due to the much smaller 
core diameter [d = 1.6 /im compared to d = 11.2 fim 
in Fig. |4]) and mode areas (ai = 0.43 /im^ compared 
to ai = 21.2 /im^ in Fig. [4|) a compression similar to 
the one shown in Fig. |4] can be done with 10 times lower 



pulse energy and power: one can compare the compressed 
pulse shown in Fig. [7] for A = 1.0 /im with the one for 
A = 7.0 /im shown in Fig. ID Both have the same soliton 
order, and achieve essentially the same level of compres- 
sion, but the pump power required is Pin = 3.1 kW when 
A = 1.0 /im and Pin = 41.8 kW when A = 7.0 /im. 

4. Conclusions 

Our recent theory for cascaded quadratic soliton com- 
pression [10,12,23] has shed new light on the role of the 
dispersion in the compression of fs near-IR pulses with 
cascaded SHG. But what can be gained by having the 
control over the GVM that microstructured fibers of- 
fer [3]? An initial study on silica microstructured fibers 
illustrated that zero GVM is not a definite advantage [14] 
because it surprisingly gave a more distorted compressed 
pulse. Moreover a very large phase mismatch was ob- 
served [2] , which required a very large effective quadratic 
nonlinear response by thermal poling of the fiber. 

Therefore we here investigated index-guiding mi- 
crostructured polymer optical fibers, because poling of 
polymers or the inclusion in the polymer matrix of nano- 
materials potentially can meet the requirement of a large 
effective quadratic nonlinear response des- A definite an- 
swer to the required nonlinearity cannot be given since 
at present the investigated polymer type (Topas) has 
yet to have its material Kerr nonlinear refractive index 
"■Kerr measured. However, based on n^j,j.j.-measurements 
of other similar polymer types like PMMA, we estimate 
that it is around 10 times that of silica. Consequently 
doff must exceed 5 pm/V and ideally be around 10 pm/V 
if efficient generation of large negative nonlinear phase 
shifts. This does pose a serious challenge for fabricating 
such an mPOF, but we believe it should be possible. 

GVM is considered a major obstacle for obtain- 
ing high-quality [12] and short duration compressed 
pulses [10]. In order to substantially reduce GVM we 
found from dispersion calculations of the transverse fiber 
modes that the relative hole size had to be quite large 
and at the same time the hole pitch had to be close to the 
wavelength of the pump light; this is consistent with the 
results obtained in silica microstructured fibers [3] . The 
problem with choosing such a design is that the SH is 
no longer single moded, and partial phase matching can 
take place to any of the SH transverse modes, thereby 
giving rise to many cascaded effects. This is clearly un- 
desirableEl On the other hand, the main goal of reducing 
GVM was to enter the stationary regime, where efficient 
and clean compression may take place. Since the fiber 
designs turned out to be in the stationary regime almost 
irrespective the size of the hole pitch or relative hole size, 
then the motivation for reducing the GVM fades. 

The main reason why the fiber designs are always in 
the stationary is the large phase mismatch. Thus, while 

*Note that the study in [3] in case of a multi-modcd SH assumed 
that phase matching to the lowest order SH mode was done using 
QPM, so a single-moded SH was not as crucial as it is here. 
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the phase mismatch gives problems in demanding a large 
doff, it has the benefit of leaving the cascaded interac- 
tion in the stationary regime, which is exactly charac- 
terized by the phase mismatch dominating over GVM 
effects [10,12]. It is worth to remark that we here present 
a fiber design where compression can occur in the sta- 
tionary regime at Ai ~ 800 nm; most bulk quadratic 
nonlinear materials are in the nonstationary regime at 
this wavelength [12]. 

The surprising conclusion is therefore that a fiber de- 
sign, where the dispersion is not too different from the 
material dispersion, gave the best compression. Numer- 
ical simulations evidenced that few-cycle pulses at a 
wavelength of 800 nm can be generated from oscillator- 
level pulse energies (a few nJ). Such pulse compression 
was demonstrated with an endlessly single-mode mPOF 
with a relative hole size of 13 = 0.4 both with a large 
pitch A = 7.0 /im and a small pitch A = 1.0 fiia. In both 
cases an effective nonlinearity of dcs ~ 10 pm/V was 
used. The only advantage of the small-pitch design was 
that much lower pump powers and pulse energies are re- 
quired, but on the other hand it is much more difficult 
to actually draw such a small-pitch mPOF. 

A natural question is whether the polymer material 
can sustain the large intensities involved when dealing 
with fs pulses; typically kWs of power are focused to 
mode areas of a few fim^ giving intensities ranging from 
10—1000 GW/cm-^. However, because the fs pulses are so 
short in time, the energy fluences in the examples shown 
in this paper were < 10 mJ/cm^, a value that unlike 
the intensity does not change much during compression. 
Such fluences should be below the threshold for writing, 
e.g., gratings in polymer materials [29], and we therefore 
expect that material damage can be avoided. 

We can therefore state that if an mPOF can be drawn 
with large enough effective nonlinearities to outbalance 
the material Kerr nonlinearity, then only a few em's of 
fiber is needed in order to generate clean few-cycle pulses 
in the near-IR with very low pulse energies. This would 
be a remarkable add-on to any fs pulsed oscillator. 
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tion equations for the electrical fields £j{z, r) read 



CiSi + KsHG,i'5'i^-ri^^2e 

+ KKcrr.l'S'l 
JC2£2 



(14a) 



£l (.Fii|fi|2+2.Fi2|f2|') 



= 



(14b) 



r,2<S'2 [^2 {^22\£2 



2Ti2\Si\')] =0 



The fields are taken scalar E(x, z, t) 



c.c, by assuming they 



are polarized along the same polarization direction 
X. The transverse field is split from the longitudinal 
propagation field by looking for solutions on the form 



ii;,(x,z,t)=£,(z,t)F,(x)e*ft- 



(15) 



where (ij are the mode propagation wave numbers and 
Fj the transverse mode profiles, and x — [x, y). 
The linear propagation operators are 



- (9 a - 
i-i = i-^ + + Di, 

oz 2 

L-2 = «Tr + ^ idl2^ + D2,cS 

oz 2 OT 
where Dj are dispersion operators up to order 



E 



m=2 
D2.c« ^b2 + 



ml (?r™ 

^-1 "12 <^ 



(16a) 
(16b) 

(17a) 



The unusual form of the SH dispersion (|17b[) is discussed 
below. The mode effective indices nes.j are related to the 
propagation constants as (3j = ricsjujj/c. The fields are 
in the frame of reference traveling with the FW group 
velocity Vg^i by the transformation t = t — z/vg^i, which 
gives the group-velocity mismatch term di2 — l/wg,i — 

l/vg^2, where Vg^, = pf\ and (3^^^ = 3^ {3, / du^^^^r 
Linear losses are included through the loss-parameter a. 
Finally, A/3 = (32 — 2/3i is the phase mismatch. 
The quadratic nonlinear coefficients are 



|/dxx(2)(x)[Fr(x)]2F2(x)| 



2cneff, 



(18) 



4ff |/^dx[J^l*(x)]2F2(x)| 



where the mode overlap area is 



A. Generalized propagation equations 

For SHG in a wave-guiding medium with both quadratic 
and cubic nonlinear material response, a general prop- 
agation equation for the electrical fields including self- 
steepening and cubic higher order nonlinear terms can 
readily be derived from Refs. [23,26]. There bulk propa- 
gation is considered, but with some slight modifications 
that will become apparent in what follows the propaga- 



dx|F,(x) 



(19) 



and x'^'' is the quadratic nonlinear tensor value along 
the polarization direction of the interacting waves, and 
deff = X*-^V2 in the reduced Kleinman notation. The 
dependence of x^^' on x implies that in an index-guiding 
mPOF the nonlinearity is only present in the polymer, as 
reminded by the subscript "P" in the integral of Eq. p^ . 
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The cubic nonlinear coefRcients are 



orr,j 



Scricgjaj ca. 



■nKc 



(20) 



where nxorrj = 3Re(x^^^)/8ncftj, and x^'^^ is the cubic 
nonUnear tensor in the polarization direction of the in- 
teracting waves. We neglect two-photon absorption so 
Ini(x(^-') = 0. The mode overlap integrals are defined as 



^dx|F,(x)nF,(. 



(21) 



In order to describe ultra-short pulses adequately, the 
usual approach adopted in silica fibers is to divide the 
material Kerr response must be divided in an electronic 
response and a vibrational (Raman) response [30] . How- 
ever, for the polymer materials used for optical fibers, 
the delayed nature of the Raman response is unknown. 
Therefore we decided to neglect the delayed vibrational 
Raman response. This can be also be justified by the 
fact that the propagation distances are very small, on 
the order of a few cm. 

Finally, we have included steepening terms through a 
self-steepening operator Sj = 1 + ^ ^ • 

The equations are valid in the slowly evolving 
wave approximation (SEWA) [25], which is a general 
spatio-temporal model with space-time focusing terms 
important for describing fs spatio-temporal optical soli- 
tons. It was recently extended to SHG by Moses and 
Wise [26] , and as a plane- wave model for SHG with com- 
peting cubic nonlinearities by Bache et al. [23]. The ad- 
vantage of the SEWA model is that it does not pose 
any constriction on the pulse bandwidth, and there- 
fore holds to the single-cycle regime. Instead, the more 
commonly used slowly varying envelope approximation 
(SVEA) only holds for Aw/tj < 1/3 (and that only when 
including steepening terms and the general Raman con- 
volution response [30]). In absence of diffraction, the dif- 
ference between the SHG SEWA model and the usual 
SVEA model is that the SH has an effective dispersion 
term (jl7bp . In the SHG SEWA model we must remember 
that one assumption made when deriving Eqs. (|14p was 
that the spectra of the fundamental and SH do not over- 
lap (substantially). This assumption allows us to sepa- 
rate the fields in two waves. We chose Aw/tjj =0.9. This 
could give some overlap between the fundamental and 
SH spectra, but we always made sure that the spectral 
components in the overlapping regions were negligible. 

We now rescale space and time (in our notation, a 
primed variable is always dimensionless) so z' = z/Ld,i, 

t' = r/Tia, where Ld.i = Tj^j/I/Sp"*! is the characteristic 
GVD length of the FW and Tin the input pulse duration. 
The fields are now normalized to the peak input electric 
field £in = £{z = 0,t = 0). Since the SH has no input 
field, we choose to rescale it to the FW input field £ia, 
so Ui = £i/£in and U2 = £il\f^£\i\ , and the equations 



become 

2n/i2 



C'2U2 + y/\W\NsHGS2U?< 



2nV22 

hi 



fll ' 

-iAfS'z' 



\U2 



2/r 
nfx 



= 



\Ui\ 



(22a) 

(22b) 

= 



where the dimensionless SHG soliton number is defined 
in Eq. ^ (an extension of the bulk soliton number from 
Refs. [9,23] to the wave-guiding case). The dimension- 
less phase mismatch is A/3' = A/3Ld,i. The cubic soliton 
number A^kcit is given by Eq. ^ and is well known from 
the nonlinear Schrodingcr equation (NLSE) in fiber op- 
tics [24] . In Eq. ((22|l the usual overlap integrals appear 

f :F,u /prfx|F,-(x)|^|ffe(x)|^ 
^^''a,a, /dx|F,(x)|2/dx|Ffe(x)|2 ^ ^ 

Finally and n = 7T.cff,i/'^cff,2 which is typically close to 
unity. The dimensionless propagation operators are 



dz' 
d_ 



Qn 



d 



m=2 



(24a) 
(24b) 
(24c) 



where we have introduced the dimensionless loss a 
aLD i and the dimensionless dispersion coefficients 



d'12 = di2 



(25) 



Finally, the steepening operators working with dr 
mensionless time are S'l = (l + is'^ 



where s' = (ti^iTjf 



) and S'^ = 
The SH effective 



dispersion (|17bp in dimensionless form is 

„ Q2 



D'2 



^2 2 9r'2 



where the dimensionless factor v ; 



ing S'^-^ = E™=o ) £^ we get [23] 



E 



Or" 
'2 



(26) 

By us- 

(27) 



The dimensionless propagation equations (|22p are the 
starting point of the analysis. The difference between the 
bulk equations we presented in Ref. [23] is that the soli- 
ton numbers and the coefficients for the SPM and XPM 
terms are modified to include the mode overlap areas, 
and that we are dealing with power and mode propa- 
gation constants instead of intensity and wave numbers. 
However, the dimensionless form is general, so the scal- 
ing laws and the critical transition points to compression 
found in Ref. [23] will still hold. 
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B. Calculation of the transverse fiber modes 
and dispersion 

The propagation equations describe the dynamics of the 
field envelope in the z propagation direction, and were 
found by describing the field as in Eq. psp . The trans- 
verse modes Fj (x) will from this analysis have to obey a 
Helmholtz-type of equation, whose solution will give the 
transverse eigenmodes Fj (x) and corresponding cigen- 
frequencies allowed by the fiber. 

We calculated the fiber modes with the MIT Photonic- 
Bands (MPB) package [31]. Each unit cell contained 
= 48^ grid points, and the super cell contained 
= 7^ unit cells. For a given (if^^ = (3-jK/2TT, the fun- 
damental mode frequency lo^^^ — loiK/c and group ve- 
locity were first calculated, followed by iterations of the 
SH until \uf^^ - 2tjf < iQ-A Material dispersion, 
parameterized by a Sellmeier equation (see App.[C]), was 
then included using a perturbative technique [32] , whose 
advantage is that many different A values can be cal- 
culated perturbatively from the MPB data (where A is 
unity). From these modified data we may then calculate 
the dispersion properties of the fiber including the effect 
of material dispersion. We should mention that after the 
perturbative technique is applied we get a modified set 
of eigenfrequencies Co^^^ , and therefore the requirement 
|(^MPB _ 2c:,MPB| ;^Q-4 lo^gQT: holds. However, the 
j-^MPB j^MPB-j (jg^^g^ gg^g -vjrere afterwards converted to di- 
mensional form, and then fitted to a regular grid. This 
ensures that the calculation of the SH dispersion was 
actually done at the proper frequency L02 — 2uji. The 
higher-order dispersion used in the numerics was calcu- 
lated with a robust polynomial fitting routine, that gave 
proper convergent results compared to the original fij- 
values when fitting up to 10 polynomial orders. 

C. Topas Sellmeier equation and losses 

We here make an accurate fit with a Sellmeier equation 
to the refractive index data from [20] made in the visible 
and near-IR. Using Mathematica we fit to the follow- 
ing single-resonance Sellmeier equation 

n2(A) = 1+B/(1- A/A^) (28) 

where A is the wavelength measured in /^m. The measure- 
ments were made from 15-75 °C and with wavelengths 
between 435.8 nm and 1014 nm [20]. The fitting param- 
eters presented here are more accurate than the ones 



T [°C] 


A [/im^] 


B 


15 


1.08199-10'^ 


2 


1.31211 


25 


1.09177-10- 


2 


1.30835 


50 


1.08080-10- 


2 


1.30139 


75 


1.09199-10" 


2 


1.29348 



Table 1 . Sellmeier equation fitting parameters for re- 
fractive index data points in Topas grade 5013 measured 
in [20] at various temperatures and with wavelengths be- 
tween 435.8 nm and 1014 nm. 




X [|am] 

Fig. 8. The loss-parameter a used in the numerics: both 
linear losses (estimated to 0.5 dB/cm [20]) and the ab- 
sorption peaks of Topas, see [20] , are modeled. 

reported in [20] (could be because in [20] a 3-parameter 
fit was used, while we use a 2-parameter fit). 

The linear losses were included in the numerical sim- 
ulations, but had only little influence since the fiber 
lengths considered here were on the order of a few em's. 
The absorption peaks of Topas around A = 1.2 /zm 
and A = 1.4 /im were modeled using a loss method: 
top-hat transmission profiles with a maximum trans- 
mission of unity were fitted to the three main spectral 
windows A S [0.29,1.21] /im, A S [1.25,1.35] /im and 
A G [1.4,1.7] /im as measured in [20] in a L = 3.2 
mm sample. The linear losses were then calculated as 
a ~ — ln(r)/L-|-Q!o, where T is the top-hat transmission 
profile, and ao are the base linear losses (estimated to 
0.5 dB/cm [20], i.e. ao = 0.115 cm"!). The loss profile 
of the simulations is shown in Fig. [51 

References 

1. T. Birks, D. Mogilevtsev, J. Knight, and P. St. J. Rus- 
sell, "Dispersion compensation using single-material 
fibers," Photonics Technology Letters, IEEE 11, 674- 
676 (1999). 

2. J. Laegsgaard, P. J. Roberts, and M. Bache, "Tailor- 
ing the dispersion properties of photonic crystal fibers," 
Opt. Quant. Electron. 39, 995-1008 (2007). 

3. M. Bache, H. Nielsen, J. Lagsgaard, and O. Bang, "Tun- 
ing quadratic nonlinear photonic crystal fibers for zero 
group-velocity mismatch," Opt. Lett. 31, 1612-1614 
(2006). 

4. G. I. Stegeman, D. J. Hagan, and L. Torner, "x'^' cas- 
cading phenomena and their applications to all-optical 
signal processing, mode-locking, pulse compression and 
solitons," Opt. Quant. Electron. 28, 1691-1740 (1996). 

5. R. DeSalvo, D. Hagan, M. Sheik-Bahae, G. Stegeman, 
E. W. Van Stryland, and H. Vanherzeele, "Self-focusing 
and self-defocusing by cascaded second-order efi'ects in 
KTP," Opt. Lett. 17, 28-30 (1992). 

6. C. R. Menyuk, R. Schiek, and L. Torner, "Solitary waves 
due to x*^^ : x'^' cascading," J. Opt. Soc. Am. B 11, 
2434-2443 (1994). 

7. X. Liu, L. Qian, and F. W. Wise, "High-energy pulse 
compression by use of negative phase shifts produced 



10 



by the cascaded x'^"* • x'^' nonlinearity," Opt. Lett. 24, 
1777-1779 (1999). 

8. S. Ashihara, J. Nishina, T. Shimura, and K. Kuroda, 
"Soliton compression of femtosecond pulses in quadratic 
media," J. Opt. Soc. Am. B 19, 2505-2510 (2002). 

9. J. Moses and F. W. Wise, "Soliton compression in 
quadratic media: high-energy few-cycle pulses with a 
frequency-doubling crystal," Opt. Lett. 31, 1881-1883 
(2006). 

10. M. Bache, O. Bang, W. Krolikowski, J. Moses, and 
F. W. Wise, "Limits to compression with cascaded 
quadratic soliton compressors," Opt. Express 16, 3273- 
3287 (2008). 

11. F. O. Ilday, K. Beckwitt, Y.-F. Chen, H. Lim, and F. W. 
Wise, "Controllable Raman-like nonlinearities from non- 
stationary, cascaded quadratic processes," J. Opt. Soc. 
Am. B 21, 376-383 (2004). 

12. M. Bache, O. Bang, J. Moses, and F. W. Wise, "Nonlo- 
cal explanation of stationary and nonstationary regimes 
in cascaded soliton pulse compression," Opt. Lett. 32, 
2490-2492 (2007). 

13. J. Moses, E. Alhammah, J. M. Eichenholz, and F. W. 
Wise, "Efficient high-energy femtosecond pulse compres- 
sion in quadratic media with flattop beams," Opt. Lett. 
32, 2469-2471 (2007). 

14. M. Bache, J. Lagsgaard, O. Bang, J. Moses, and F. W. 
Wise, "Soliton compression to ultra-short pulses using 
cascaded quadratic nonlinearities in silica photonic crys- 
tal fibers," Proc. SPIE 6588, 65880P (2007). 

15. D. Faccio, A. Busacca, W. Belardi, V. Pruneri, P. Kazan- 
sky, T. Monro, D. Richardson, B. Grappe, M. Cooper, 
and C. Pannell, "Demonstration of thermal poling in 
holey fibres," Electron. Lett. 37, 107-108 (2001). 

16. P. G. Kazansky, private communication. 

17. P. G. Kazansky and V. Pruneri, "Electric-field poling of 
quasi-phase-matched optical fibers," J. Opt. Soc. Am. B 
14, 3170-3179 (1997). 

18. C. Chang, C. Chen, C. C. Chou, W. J. Kuo, and J. J. 
Jeng, "Polymers for electro-optical modulation," Jour- 
nal of Macromolecular Science: Polymer Reviews 45, 
125-170 (2005). 

19. G. Y. Guo and J. C. Lin, "Second-harmonic generation 
and linear electro-optical coefficients of BN nanotubes," 
Phys. Rev. B 72, 075416 (2005). 

20. G. Khanarian and H. Celanese, "Optical properties 
of cyclic olefin copolymers," Opt. Eng. 40, 1024-1029 
(2001). 

21. L. F. MoUenauer, R. H. Stolen, and J. P. Gordon, "Ex- 
perimental observation of picosecond pulse narrowing 
and solitons in optical fibers," Phys. Rev. Lett. 45, 
1095-1098 (1980). 

22. W. Krolikowski, O. Bang, N. Nikolov, D. Neshev, J. 
Wyller, J. Rasmussen, and D. Edmundson, "Modula- 
tional instability, solitons and beam propagation in spa- 
tially nonlocal nonlinear media," ,J. Opt. B: Quantum 
Semiclass. Opt. 6, s288 (2004). 

23. M. Bache, J. Moses, and F. W. Wise, "Scaling laws for 
soliton pulse compression by cascaded quadratic nonlin- 
earities," J. Opt. Soc. Am. B 24, 2752-2762 (2007). 

24. G. P. Agrawal, Nonlinear fiber optics , 3 ed. (Academic 
Press, London, 2001). 

25. T. Brabec and F. Krausz, "Nonlinear Optical Pulse 



Propagation in the Single-Cycle Regime," Phys. Rev. 
Lett. 78, 3282-3285 (1997). 

26. J. Moses and F. W. Wise, "Controllable Self-Steepening 
of Ultrashort Pulses in Quadratic Nonlinear Me- 
dia," Phys. Rev. Lett. 97, 073903 (2006), see also 
arXiv:physics/0604170. 

27. B. T. Kuhlmey, R. C. McPhedran, and C. M. de Sterke, 
"Modal cutoff in microstructured optical fibers," Opt. 
Lett. 27, 1684-1686 (2002). 

28. T. A. Birks, J. C. Knight, and P. S. RusseU, "Endlessly 
single-mode photonic crystal fiber," Opt. Lett. 22, 961- 
963 (1997). 

29. A. Baum, P. J. Scully, W. Perrie, D. Jones, R. Is- 
sac, and D. A. Jaroszynski, "Pulse-duration dependency 
of femtosecond laser refractive index modification in 
poly(methyl methacrylate) ," Opt. Lett. 33, 651-653 
(2008). 

30. K. J. Blow and D. Wood, "Theoretical description of 
transient stimulated Raman scattering in optical fibers," 
IEEE J. Quantum Electron. 25, 2665-2673 (1989). 

31. S. Johnson and J. Joannopoulos, "Block-iterative 
frequency-domain methods for Maxwell's equations in 
a planewave basis," Opt. Express 8, 173-190 (2001). 

32. J. Lffigsgaard, A. Bjarklev, and S. Libori, "Chromatic 
dispersion in photonic crystal fibers: fast and accurate 
scheme for calculation," J. Opt. Soc. Am. B 20, 443-448 
(2003). 



11 



